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ABSTRACT 

We demonstrate that high abundances of water vapor could have existed in extremely low metallicity (10^^ 
solar) partially shielded gas, during the epoch of first metal enrichment of the interstellar medium of galaxies 
at high redshifts. 
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1. INTRODUCTION 

Water is an essential ingredient for life as we know it 
jKastingl l2010t) . In the interstellar medium (ISM) of the 
Milky-Way and also in external galaxies, water has been ob¬ 
served in the gas phase and as grain surface ice in a wide vari¬ 
ety of environments. These environments include diffuse and 
dense molecular clouds, photon-dominated regions (PDRs), 
shocked gas, protostellar envelopes, and disks (see review by 
Ivan Dishoeck et aD2013l) . 

In diffuse and translucent clouds, H2O is formed 
mainly in gas phase react ions via ion-molecule sequences 
(iHerbst & Klemnere^l 19731) . The ion-molecule reaction net¬ 
work is driven by cosmic-ray or X-ray ionization of H and H2, 
that leads to the formation of H+ and H^ ions. These inter¬ 
act with atomic oxygen and form OH+. A series of abstrac¬ 
tions then lead to the formation of H3O+, which makes OH 
and H2O through dissociative recombination. This formation 
mechanism is generally not very efficient, and only a small 
fraction of the oxygen is converted into wate r, the rest remains 
in ato mic form, or freezes out as water ice (iHollenbach et al.l 

[2^. __ 

ISonnentrucker et ^ (120101) showed that the abundance of 
water vapor within diffuse clouds in the Galaxy is remarkably 
constant, with xhjO 10^^, that is ~ 0.1% of the available 
oxygen. Here ahjO is the H2O number density relative to the 
total hydrogen nuclei number density. Towards the Galactic 
center this value can be enhanced by up to a factor of ~ 3 
dMonie et al.ll201 ll ISonnentrucker et £0120131) . 

At temperatures > 300 K, H2O may form directly via the 
neutral-neutral reactions, O H- H2 —OH + H, followed by OH 
H- H2 —!> H2O + H. This formation route is particularly impor¬ 
tant in shocks, where the gas heats up to high t emperatures, 
and can drive most of the oxyg en into H2O (iDraine et al.l 
fl^lKaufman & Neufeldlflom 

Temperatures of a few hundreds K are also expected in very 
low metallicity gas environments, with elemental oxygen and 
carbon abundances of < 10~^ solar (iBromm & LoebI 120031: 
lOmukai et al.l 120051: iGlover & ClarkI 120141) . associated with 
the epochs of the first enrichment of the ISM with heavy el- 
ements, in the first gene ration of galaxies at high redshifts 
(iLoeb & Furlanettoll2013l) . At such low metallicities, cooling 
by fine structure transitions of metal species such as the [CII] 
158 jim line, and by rotational transitions of heavy molecules 
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such as CO, becomes inefficient and the gas remains warm. 

Could the enhanced rate of H2O formation via the neutral- 
neutral sequence in such warm gas, compensate for the low 
o xygen abundance at l ow metallicities? 

lOmukai et akl (120051) studied the thermal and chemical evo¬ 
lution of collapsing gas clumps at low metallicities. They 
found that for models with gas metallicities of 10^^ — 10^^ 
solar, ahjO reach 10^^, but only if the density, n, of the 
gas approaches 10^ cm^^. Photodissociation of molecules by 
far-ultraviolet (FUV) radiation was not included in their study. 
While at solar metallicity dust-grains shield the interior of gas 
clouds from the FUV radiation, at low metallicities photodis¬ 
sociation by FUV becomes a major removal process for H2O. 
H2O photodissociation produces OH, which is then itself pho- 
to dissociated into atomic o xygen. 

[Hollenbach et alJ (120121) developed a theoretical model to 
study the abundances of various molecules, including H2O, in 
PDRs. Their model included many important physical pro¬ 
cesses, such as freezeout of gas species, grain surface chem¬ 
istry, and also photodissociation by FUV photon s. However 
they focused on solar metallicity. Intriguingly, iBavet et ^ 
(I2OO9I) report a water abundance close to 10^^ in the optically 
thick core of their single PDR model for a low metallicity of 
10^^ (with n = 10^ cm^^). However, Bayet et al. did not 
investigate the effects of temperature and UV intensity varia¬ 
tions on the water abundance in the low metallicity regime. 

Recently, a comprehensive study of molecular abundances 
for the bulk ISM gas as functio ns o f the metallicity, were stud- 
ied by iPenteado et alJ (12014l) and iBialv & SternbergI (12014 
hereafter BS15), these models however, focused on the “low 
temperature” ion-molecule formation chemistry. 

In this Letter we present results for the H2O abundance in 
low metallicity gas environments, for varying temperatures, 
FUV intensities and gas densities. We demonstrate the impor¬ 
tance of the onset of the neutral-neutral formation sequence, 
and explore the role of the FUV field. We find that for tem¬ 
peratures T in the range 250 — 350 K, H2O may be abundant, 
comparable to or higher than that found in diffuse Galactic 
clouds, provided that the FUV intensity to density ratio is 
smaller than a critical value. In § |2] we discuss our physi¬ 
cal and chemical model. We present our results in §[3] Finally 
we summarize and discuss our conclusions in §|4] 

2. MODEL INGREDIENTS 

We calculate the abundance of gas-phase H2O for low 
metallicity gas parcels, that are exposed to external FUV ra¬ 
diation and cosmic-ray and/or X-ray fluxes. Given our chemi¬ 
cal network we solve the steady state rate equations using our 
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dedicated Newton-based solver, and obtain xu^o as function 
of T and the FUV intensity to density ratio. 

We adopt a 10^ K diluted blackbody spectrum, representa¬ 
tive of radiation produced by massive Pop-Ill stars. The pho¬ 
ton density in the 6-13.6 eV interval, is Uy = nyflluy, where 
Hy.o = 6.5 X 10^^ photons cm~^ is the photon density in the 
interstellar radiation field (ISRF, lDrainel201 Ih . and /(jv is the 
“normalized intensity”. Thus /uv = 1 corresponds to the FUV 
intensity in the Draine ISRF. 

Cosmic-ray and/or X-ray ionization drive the ion-molecule 
chemical network. We assume an ionization rate per hydro- 

gen nucleon C (s~ ^). _ In the Galaxy, iDalgarnol (l2006h and 

llndriolo & McCalll (120121) showed that lies within the rel¬ 
atively narrow range 10^'^ — 10^^® s^^ We therefore intro¬ 
duce the “normalized ionization rate” = (i^/lO^*® s^*). 
The ionization rate and the FUV intensity are likely corre¬ 
lated, as both should scale with the formation rate of massive 
OB stars. We thus set C-i6 = ^uv as our fiducial case but also 
consider the cases i^_i6 = lO^'/yv and = 10/uv- 

Dust shielding against the FUV radiation becomes inef¬ 
fective at low metallicities. However, self absorption in the 
H 2 lines may significantly attenua te the destructive Lyman 
Werner (11.2-13.6 eV ) radiation (iDraine & Bertoldil 119961: 
ISternberg et al.l l?014l) and high abundances of H 2 may be 
maintained even at low metallicity (BS15). In the models pre¬ 
sented here we assume an H 2 shielding column of at least 
5 X 10^' cm^^. (For such conditions CO is also shielded by 
the H 2 .) The LW flux is then attenuated by a self-shielding 
factor of /shield ~ 10^^ and the H 2 photodissociation rate is 
only 5.8 x 10^^^/uv With this assumption H 2 photodis¬ 
sociation by LW photons is negligible compared to cosmic- 
ray and/or X-ray ionization as long as /yv < 85i^_i6. 

However, even when the Lyman Werner band is fully 
blocked, OH and H 2 O are photodissociated because their en¬ 
ergy thresholds for photodissociation are 6.4 and 6 eV, respec¬ 
tively. For the low metallicities that we consider, photodisso¬ 
ciation is generally the dominant removal mechanism for H 2 O 
and OH. We adopt the calculated OH and H 2 O photodissoci¬ 
ation rates calculated by BS15. 

We assume thermal and chemical steady states. In the 
Milky Way, the bulk of the ISM gas is considered to be at 
approximate thermal equilibrium, set by cooling and heating 
processes . We discuss the relevant chemical and thermal time- 
scales in (12.21 

Given the above mentioned assumptions, the steady state 
solutions for the species abundances depend on only two pa¬ 
rameters, the temperature T and the intensity to density ratio 
^uv/« 4 - Here 114 = cm^^) is the total number den¬ 

sity of hydrogen nuclei normalized to typical molecular cloud 
densities. T and /uv/n 4 form our basic parameter space in 
our study. 


2.1. Chemical Model 

We consider a chemical network of 503 two-body reac¬ 
tions, between 56 atomic, molecular and ionic species of H, 
He, C, O, S, and Si. We assume cosmological elemental he¬ 
lium abundance of 0.1 relative to hydrogen (b y number). For 
the me tal elemental abundances we adopt the lAsplund et al.l 
(l2009l) photospheric solar abundances, multiplied by a metal¬ 
licity factor Z' (i.e., Z' = 1 is solar metallicity). In our fiducial 
model we assume Z' = 10^^, but we also explore cases with 
Z' = 10^^ and Z' = 10^^. Since our focus here is on the very 


low metallicity regime, where dust grains play a lesser role, 
we neglect any depletion on dust grains, and dust-grain chem¬ 
istry (except for H 2 , as discussed further below). Direct and 
induced ionizations and dissociations by the cosmic-ray / X- 
ray field (oc Q are included. For the gas phase reactions, we 
adopt the rate coeffici ents given by the UMIST 2012 database 
(iMcElrov et alJl2013h . 

The formation of heavy molecules relies on molecular hy¬ 
drogen. We consider two scenarios for H 2 formation: (a) pure 
gas-phase formation; and (b) gas-phase formation plus forma¬ 
tion on dust grains. In gas-phase, radiative-attachment 

H -f e ^ H- -b V (1) 

followed by associative-detachment 

H + H ^ H 2 + e (2) 

is the dominant H 2 formation route. 

H 2 formation on the surface of dust grains is considered 
to be the dominant formation channel in the Milky way. We 
adopt a standard rate coefficient 

R ~ 3 X 10^'^ Z! cm^ s^' (3) 

(iHollenbach et alJll971l : Llurall 19741 : ICazaux & Tielensll2002l) . 

where T 2 = ij /lOO K). In this expression we assume that 
the dust-to-gas ratio is linearly proportional to the metallicity 
Z'. Thus, in scenario (b) H 2 formation on dust grains domi¬ 
nates even for Z! = 10^^. Scenario (a) is the limit where the 
gas-phase channel dominates, as appropriate for dust-free en¬ 
vironments or for superlinear dependence of the dust-to-gas 
ratio on Z' (see also BS15). 


2.2. Time scales 

The timescale for the system to achieve chemical steady 
state is dictated by the relatively long H 2 formation timescale. 
In the gas phase (scenario [a]) it is 

tR, = 8 X 10^ «4yr - (4) 

k 2 nxQ 10 ^ z. 


where Xg = 2.4 x 10 is the electron frac¬ 

tion as set by ionization recombination equilibrium, and ki ~ 

l. 5 X cm^ s^' is the rate coefficient for reaction 

m. 

For formation on dust grains ( 0 = Z'), the timescale is gener¬ 
ally shorter, with 


tR 2 


1 

Rn 


10“ % 


- 1/2 



yr. 


(5) 


Gas clouds with lifetimes t ^ t] 4 ^ will reach chemical steady 
state. 

The relevant timescale for thermal equilibrium is the cool¬ 
ing timescale. For low metallicity gas with Z' = 10^^, 
the cooling proceeds mainly via H 2 rotational excitations 
(iGlover & ClarkI 12014ft . If the cooling rate per H 2 molecule 
(in erg s^^) is W(n, T), then the cooling timescale is given by 

_ ksT 

W{n,T) ■ ^ ’ 

Here ks is the Boltzmann constant. For n = 10"* cm~^, and 
T = 300 K, W ~5 X 10^^^(xh2/0.1) ergs^^ (iLe Bourlot et alJ 
119991) . and the cooling time is very short w 2 x 10^(0 .I/xhj) 
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Fig. 1.— The fractional H 2 O abundance xh^o as a function of T and /uv/ti 4 > for Z' = 10 ^ and f-ie/Ajv = 1> assuming pure gas-phase chemistry (scenario 
[a] - left panel), and including H 2 formation on dust grains (scenario [b] - right panels). In both panels, the solid line indicates the 10“® contour, which is 
a characteristic value for the H 2 O gas phase abundance in diffuse clouds in the Milky-Way. At high temperatures (or low /uv/w 4 values) the neutral-neutral 
reactions become effective and xjjjO rises. 


yr. For densities much smaller than the critical density for 
H 2 cooling, W n and fcooi !/”■ '^he opposite limit, W 
saturates and fcooi becomes independent of density. We see 
that generally fcooi < fuj- 

Because the free fall time 

is generally much shorter than chemical steady state may 
be achieved only in stable, non-collapsing clouds, with life¬ 
times ^tff. Obviously both fHj and fcooi must be also shorter 
than the Hubble time at the redshift of interest. 

3. RESULTS 

Next we present and discuss our results for the steady state, 
gas-phase H 2 O fraction XH 2 O = «H 2 o/«^ as function of tem¬ 
perature T, and the FUV intensity to density ratio /yv /fi 4 - 

3.1. XH 2 O as a function of T andluy/n^ 

Figure [U shows logi n(xH 2 o) contours for the two scenar¬ 
ios described in § 12.ll In one H 2 forms in pure gas-phase 
(scenario [a] - left panel), and in the other H 2 forms also on 
dust-grains (scenario [b] - right panel). Our fiducial parame¬ 
ters are Z' — 10^^ and i^_i6 = /uv- At the upper-left region 
of the parameter space, XH 2 O is generally low < 10^^. In this 
regime, H 2 O forms through the ion-molecule sequence, that 
is operative at low temperatures. In the lower right corner, the 
neutral-neutral reactions become effective and XH 2 O rises. 

In both panels, the solid line highlights the XH 2 O = 10 ^ 
contour, which resembles the H 2 O gas phase abundance 
in diffuse and translucent Galactic clouds. This line also 
delineates the borderline between the regimes where H 2 O 
forms via the “cold” ion-molecule sequence, and the “warm” 
neutral-neutral sequence. The temperature range at which 
the neutral-neutral sequence kicks-in is relatively narrow ^ 
250 — 350 K, because the neutral-neutral reactions are limited 
by energy barriers that introduce an exponential dependence 
on temperature. 


The dependence on /yy /m 4 is introduced because the FUV 
photons photodissociate OH and H 2 O molecules and there¬ 
fore increase the removal rate of H 2 O and at the same time 
suppress formation via the OH + H 2 reaction. This gives rise 
to a critical value for /yy /n 4 below which H 2 O may become 
abundant. 

For pure gas-phase H 2 formation (scenario [a] - left panel), 
the gas remains predominantly atomic, and H 2 O formation is 
less efficient. In this case XH 2 O ^ 10^^ only when /uv/n 4 is 
smaller than a critical value of 

(/uv/«4)i?|t = 2 X 10-2 . (8) 

However, when H 2 formation on dust is included (scenario 
[b] - right panel), the hydrogen becomes fully molecular, and 
H 2 O formation is then more efficient. In this case xhjO iriay 
reach IQ-^ for/uv/n 4 smaller than 

(/uv/«4)iril = 3 X 10-' , (9) 

an order of magnitude larger than for the pure-gas phase for¬ 
mation scenario. 

3.2. Variations in Z'and 

In Figure I 2 ] we investigate the effects of variations in the 
value of Z' and the normalization C-ie/^uv- The Figure 
shows the XH 2 O = 10-^ contours for scenarios (a) (left panels) 
and (b) (right panels). As discussed above, H 2 O is generally 
more abundant in scenario (b) because the hydrogen is fully 
molecular in this case, and therefore the 10-^ contours are 
located at higher I\jy jn^^. values in both right panels. 

The upper panels show the effect of variations in the metal- 
licity value Z', for our fiducial normalization i^_i6 = /(jy. In 
both panels, the oxygen abundance rises, and XH 2 O increases 
with increasing Z'. Thus at higher Z', the IQ-^ contours shift 
to lower T and higher /yy /n 4 and vice versa. An exception is 
the behavior of the Z' = IQ-^ curve, for which the metallic- 
ity is already high enough so that reactions with metal species 
dominate H 2 O removal for /uv /«4 ^ IQ-^. The increase in 
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Fig. 2.— The xhjO = 10 * contour, for variations in 7! (upper panels) and in 
panels), and including H 2 formation on dust grains (scenario [b] - right panels). 

metallicity then results in a decrease of the H 2 O abundance, 
and the 10^^ contour shifts to the right. For I\]\ jiiA, > 10^^ 
removal by FUV dominates and the behavior is similar to that 
in the Z' = 10^^ and Z' = 10^"^ cases. 

The lower panels show the effects of variations in the ion¬ 
ization rate normalization C-ie/Zuy for our fiducial metallic¬ 
ity value of Z' = 10^^. First we consider the pure gas phase 
formation case (scenario [a] - lower left panel). For the two 
cases i^_i6//uv = 1 and 10^*, the H 2 O removal is dominated 
by FUV photodissociation and therefore is independent of 1 ^. 
As shown by BS15, the H 2 O formation rate is also indepen¬ 
dent of when the H 2 forms in the gas-phase. Therefore 
JCH 2 O is essentially independent of 1 ^ and the contours over¬ 
lap. For the high ionization rate C-ie/fuv = 10, the proton 
abundance becomes high, and H 2 O reactions with H+ domi¬ 
nate H 2 O removal. In this limit xh,o decreases with 1 ^ and 
the 10^^ contour moves down. 

When H 2 forms on dust (scenario [b] - lower right panel), 
the H 2 O formation rate via the ion-molecule sequence is pro¬ 
portional to the H+ and abundances, which rise with 
Since the gas is molecular, the proton fraction is low and the 
removal is always dominated by FUV photodissociations (in¬ 
dependent of Q. Therefore, in this case XH 2 O increases with 
C-ie/fuv and the curves shift up and to the left, toward lower 
T and higher /yv jn^- 

4. SUMMARY AND DISCUSSION 

In this Letter we have demonstrated that the H 2 O gas phase 
abundance may remain high even at very low metallicities of 
Z' ~ 10^^. The onset of the efficient neutral-neutral formation 
sequence at T ~ 300 K, may compensate for the low metal¬ 
licity, and form H 2 O in abundance similar to that found in 
diffuse clouds in the Milky Way. 

We have considered two scenarios for H 2 formation, repre¬ 
senting two limiting cases, one in which H 2 is formed in pure 
gas phase (scenario [a]), and one in which H 2 forms both in 
gas-phase and on dust grains, assuming that the dust abun¬ 
dance scales linearly with Z' (scenario [b]). Recent studies 


C-16/4jv (lower panels), assuming pure gas-phase chemistry (scenario [a] - left 

byjGalametzet^ (1201 Ih . iHerrera-Camus et'H] (1201 2h and 
iFisher et al.l (12014 1. suggest that the dust abundance might 
decrease faster then linearly with decreasing Z'. As shown 
by BS15, forZ' = 10^^ and dust abundance that scales as Z'^ 
with p >2, H 2 formation is dominated by the gas phase for¬ 
mation channel. Therefore our scenario (a) is also applicable 
for models in which dust grains are present, with an abun¬ 
dance that scales superlinearly with Z'. For both scenarios 
(a) and (b) we have found that the neutral-neutral formation 
channel yields XH 2 O ^ 10^^, provided that /uv /«4 is smaller 
than a critical value. For the first scenario we have found that 
this critical value is (/uv/« 4 )crit = 2 x 10^. For the second 
scenario (/uv/« 4 )crit = 3 x 10^*. 

In our analysis we have assumed that the system had 
reached a chemical steady state. For initially atomic (or 
ionized) gas, this assumption offers the best conditions for 
the formation of molecules. Flowever, chemical steady state 
might not always be achieved within a cloud lifetime or even 
within the Hubble time. The timescale to achieve chemical 
steady state (from an initially dissociated state) is dictated 
by the H 2 formation process, and is generally long at low 
metallicities. For Z' = 10^^ and our fiducial parameters, the 
timescal es for both scena rios are of order of a few 10^ years 
(see e.g. iBell et al.ll20d^ . and are comparable to the age of 
the Universe at redshifts of ^ 10. The generically high water 
abundances we find for warm conditions and low metallicities 
will be maintained in dynamically evolving systems so long 
as they remain H 2 -shielded. 

Our results might have interesting implications for the ques- 
tion of how early could have life originated in the Universe 
(lLoebll20T^ . Our study addresses the first step of H 2 O for¬ 
mation in early enriched, molecular gas clouds. If such a 
cloud is to collapse and form a protoplanetary disk, some of 
the H 2 O may make its way to the surfaces of forming plan¬ 
ets (Ivan Dishoeck et aH 12014ft . However, the question of to 
what extent the H 2 O molecules that were formed in the ini¬ 
tial molecular clouds, are preserved all the way through the 
process of planet formation, is beyond the scope of this paper. 
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